Generalized calculation of magnetic coupling constants for Mott-Hubbard insulators: 

Application to ferromagnetic Cr compounds 
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Using a Rayleigh-Schrodinger perturbation expansion of multi-band Hubbard models, we present 
analytic expressions for the super-exchange coupling constants between magnetic transition metal 
ions of arbitrary separation in Mott-Hubbard insulators. The only restrictions are i) all ligand ions 
are closed shell anions and ii) all contributing interaction paths are of equal length. For short paths, 
our results essentially confirm the Goodenough-Kanamori- Anderson rules, yet in general there does 
not exist any simple rule to predict the sign of the magnetic coupling constants. The most favorable 
situation for ferromagnetic coupling is found for ions with less than half filled d shells, the (relative) 
tendency to ferromagnetic coupling increases with increasing path length. As an application, the 
magnetic interactions of the Cr compounds Rb2CrCl4, CrCl3, CrBr3 and (M3 are investigated, all 
of which except CrCl3 are ferromagnets. 

71.10.Fd, 71. 20. Be, 75.10.Jm, 75.30.Et, 75.50.Dd 

I. INTRODUCTION 

The understanding of the magnetic properties of the 
nonmetallic 3d-transition metal compounds is not satis- 
fying though they have been investigated intensively ex- 
perimentally and theoretically. In particular there exists 
only a qualitative understanding of magnitude and sign 
of the magnetic coupling on the basis of the so-called 
Goodenough-Kanamori-Anderson (GKA) rules. In this 
paper we present a more systematic way of calculating 
the magnetic coupling constants of these materials based 
on the concept of super-exchange. We also test our meth- 
ods by applying them to a sub-class of the rare class of 
ferromagnetic 3d insulators, i.e. to Cr compounds, where 
the GKA rules apparently do not apply. 

In nonmetallic, magnetic 3d-transition metal oxides 
and halides the spins are localized in the incompletely 
filled 3d shells of the transition metal sites, hence-forward 
called magnetic sites. This is a consequence of large, cor- 
relation effects. They are Mott-Hubbard insulatorsEl and 
are thus described by a Hubbard-type model of the form 
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Here ef is the energy of a d orbital a at the magnetic 
site i and ej that of an orbital r at the ligand site 
I. tfj 7 ' represents the hopping from the orbital on to 
Ti, tJiT that from 77 to r' v . U a p and I a p are Coulomb 
and exchange interactions at the magnetic sites, re- 
spectively. In the atomic limit, I a @ < leads to the 
first Hund's rule (maximum spin) and U aa > U a p to 
the second Hund's rule (maximum orbital momentum). 

I a p = I a p 4<r(*)4cr' {i) c ao'(i)cpv(j) will be called Hund 
interaction in the following. Its effect is an exchange of 
spins between orbitals a and /?. The application of / to 
symmetric (anti-symmetric) spin-states therefore has the 
effect of a multiplication with +(— )I a p. In Eq. (Q) we 
neglect the Coulomb interaction at the ligand ions and 
between different sites. We will also use the approxima- 
tion U aa = U a /3 = U and I a p = Jh- These restrictions 
could be loosened but we do not expect qualitatively dif- 
ferent results. We note that the atomic part of the Hamil- 
tonian (|l|) has been used very satisfactorily in ligand field 
theory for, e,g. the interpretation of optical data of such 
compounds .□ 

Anderson has introduced the concept of super- 
exchange where the hopping of spins between the dif- 
ferent sites leads to the coupling of n the total spins Si 
and S2 of the magnetic sites 1 and 2B As is known from 
experimental resultsu'El the interaction is described very 
well by the Heisenberg Hamiltonian 



H c ff — J12 Si • S 2 



(2) 



(T,a' 



The interpretation of the interaction due to kinetic ex- 
change (different from Coulomb exchange) is reflected in 
the formulation of the GKA rules.El In a simplified version 
they read 

A) If the hopping of a spin from a singly occu- 
pied £1 of the ion Mi to the singly occupied 
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orbital (2 of the ion M2 is possible, the inter- 
action between the spins in £1 and (2 is anti- 
ferromagnetic. 

B) When a spin can hop from the singly occu- 
pied orbital £i to an unoccupied orbital i] of 
M2, a ferromagnetic contribution to the inter- 
action results, which is weaker than the anti- 
ferromagnetic one from rule A. 

With the use of these rules the sign of the coupling con- 
stant can be predicted in many cases. Especially, the 
GKA rules explain why most of the non-metallic 3d- 
transition-metal compounds do not show a macroscopic 
magnetic moment. This is related to the fact that orbital 
ordering is usually needed for a ferromagnetic interaction 
to occur. One prominent example are the CuF2-planes 
of K2CUF4, where the Jahn- Teller effect causes orbital 
order, of alternating z 2 — y 2 and z 2 — x 2 hole orbitals, 
which completely suppresses the usually dominant anti- 
ferromagnetic interaction. 

However, in many cases the GKA rules are not helpful 
to decide which one of two competing interactions dom- 
inates. Such a case is Rb2CrCi4, where, again due to 
the Jahn- Teller effect, an analogous orbital ordering as 
in K2CUF4 occurs. Nevertheless, the nearest neighbor, 
singly occupied 3x 2 — r 2 and 3y 2 — r 2 orbitals exhibit 
a strong anti-ferromagnetic coupling via p a ligand or- 
bitals, completely absent in the case of K2CUF4. The 
fact, therefore, that Ri^CrCLt is also a ferromagnet, is 
by no means a simple consequence of the GKA rules. 
Similarly, there is no obvious GKA argument for the fact 
that the isostructural tri-halides CrBr3 and Cri3 are fer- 
romagnets. 

In our paper we present a more detailed theory of 
super- exchange based on the Hamiltonian ([!]) which al- 
lows the calculation of both anti-fcrromagnetic and fer- 
romagnetic coupling constants for hopping paths via ar- 
bitrarily many ligands. 

We test our theory by applying it to the ferromagnetic 
insulators mentioned above and show that a reasonable 
choice of parameters in (|l]) will lead to magnetic coupling 
constants which agree with experiment both in magni- 
tude and sign. 

In section II we use perturbation theory to determine 
the coupling constants for pairs of d-ions which are linked 
via an arbitrary number of ligand ions. First we present 
the interaction of spin-i pairs, where the spins are ex- 
changed along a chain of ligands each possessing one un- 
occupied orbital. We investigate both the direct super- 
exchange between two singly occupied d orbitals £ and 
the two indirect processes where spin- free (empty or dou- 
bly occupied) d orbitals at the magnetic sites are in- 
volved. We then generalize these results to the case 
of spin exchange via several paths of the same order - 
here we distinguish between equivalent and inequivalent 
paths. Finally we give an expression for the coupling of 
two d-ions with arbitrary (integral) occupation which are 
coupled via P different paths with N unoccupied ligands 



each. The application of these results to Rd^CrCU and 
the Cr tri-halides is discussed in section III. 

Throughout the paper we will adopt the hole picture 
since the ligands usually are anions with a filled p shell. 

II. SUPER-EXCHANGE 

In perturbation theory the coupling between the total 
spins Si and S2 of two magnetic ions results from the ex- 
change of the spins in the singly occupied orbitals. Since 
we do not consider spin-dependent interactions like spin- 
orbit coupling in (|lj), the spin Hamiltonian is isotropic 
and in the lowest order we get 

H cff = J CiC2 S Ci ' S C 2 - ( 3 ) 
C1C2 

Here, £i and C2 denote the half-filled orbitals at the sites 1 
and 2. Assuming a strong Hund interaction, the ground- 
state has maximum spin and the replacement 



appears to be justified. Eq. (|3|) can then be written as 

"•.°(g^) 4<s s : > s < 2 S2> <•) 

j Si ■ S 2 

4 < Si >< S 2 >■ 

We begin with the coupling of spin-^ pairs. The gener- 
alization to pairs with larger total spins is possible with 
the help of Eq. (g). 

We distinguish three cases : 
A: Direct interaction of two singly occupied orbitals £1 
and C2 without the involvement of another d orbital. In 
general, the resulting coupling will be anti-ferromagnetic, 
the coupling constant will therefore be called J^r 2 - 
FE: Indirect interaction of £1 and £2 via one ore more 
empty d orbitals rj involving the Hund interaction. The 
coupling constant is always negative and will be called 

jFE 

FD: Indirect coupling of £1 and £2 via one ore more dou- 
bly occupied d orbitals involving the Hund interaction. 
The coupling constant is negative, in general, and will be 
called j££ a . 

A. super-exchange via one ligand 

We first consider the direct coupling of two equivalent 
magnetic ions 1,2 with singly occupied orbitals £1, C2 
via one ligand with an empty hole orbital, see Fig. |l| a). 
Fourth order perturbation theory is the lowest contribut- 
ing order. There exist six possibilities of exchanging the 
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spins which differ in the sequence of spin exchange. In 
the following these sequences will be called channels. For 
each channel three intermediate states m exist. Two of 
them contain only singly occupied sites (one spin at the 
ligand and the other at one of the magnetic ions). In the 
third state the two spins occupy the same ion. The exci- 
tation energy of the singly occupied states is the charge 
transfer energy (A^ 2 ) when the magnetic ion 1 (2) 

is singly occupied. We have = e T — , when the lig- 
and is empty, which is the usual case in the hole picture. 
Note, that when transforming to the electron picture, the 
numerical value of A^ remains unchanged, but the or- 
bital energies e have to be renormalized. If one of the 
magnetic ions is doubly occupied, the excitation energy 
is U ± (Afj — Af 2 ), if it is the ligand we have A^ + A^ 2 . 
The matrix elements between the states are the hopping 
integrals given by ([!]) and each channel yields the same 
factor (t (l ) 2 (i c J 2 . Note that the energy denominators 
(Eq — E m )~ l are negative, but an additional minus sign 
occurs because of the exchange of the spins. 
The total contribution of the six channels is 
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A 2 (l U-UbC, A l U + u; ClC2 
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(6) 
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C2 \A Cl 



1 



K2 



t 2 t 2 

Hi H2 



with w^fj = Af 2 — Af t . H x (H 2 ) is the hopping integral 
between orbitals (1 (£2) and the ligand. 

Assuming A^ ~ A,c_ = A^ in Eq. (^), we obtain the 
well known expression!!! 



.v.i 

C1C2 



A, 2 



1 
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With the replacements 
1 

U + 

ki t<2 

A L 

this leads to 



= Ai 



C1C2 



(7) 



(8a) 



(8b) 



configuration. Then an inter-orbital on-site exchange can 
take place, leading to a negative coupling constant. 

The intermediate high spin state can occur in two 
ways. FE: If the spin-free orbital at the magnetic ion, 
where the direct exchange takes place, is empty, the spin 
from the half-filled orbital of the other magnetic ion can 
hop into it. FD: If the spin-free orbital is doubly occu- 
pied, one of the spins can hop to the empty ligand, ener- 
getically preferred is the one anti-parallel to the ground 
state spin. Both cases lead to a polarization of the ligand 
which is parallel to the spin of the magnetic ion in case 
FE and anti-parallel to it in the case FD. We will begin 
with the investigation of FE. 

We assume a system of two magnetic ions Mi and M 2 
with only one singly occupied orbital £ and one empty 
orbital 77 each. Also, the intermediate ligand is empty, 
see Fig. @b). 

There exist only two possibilities for exchanging the 
spins. For each, the spin from the half filled orbital ( 
of one magnetic ion has to migrate to the empty or- 
bital 77 of the other ion, where the on-site exchange of 
the spins takes place, leading to the preference of the 
ferromagnetic state. Simplifying A^ = A^ 2 



A^ and 



A^ the coupling constant is given by 



FE,1 



Jh 



(U (A c - A„))2 (A c )2 
[(*Ci t m) 2 + (k2 t Vl ) 2 ] ■ 



(10) 



Hi (tvi ) represents the hopping-integral from the orbital 
Ct (Vi) to the ligand. We note that J FE is proportional 
to Jh- 



Using the simplification A^ 



jFE.l 

H1C2 



1 Jh 



UA L 2 U 



[(Hi h 



Al we have 



(*C2 tmY 



11) 



With the substitutions 



(ft, 



Hi Si 



Vi, 



(12a) 



(12b) 



J ciC2 = Al (kt( 



(9) 



i represents the effective hopping from one to the other 
magnetic ion via the ligand. J A is always non-negative 
and leads, because of the Pauli principle, to an anti- 
parallel ordering of Si and S2- 

For the indirect coupling of the spins at the magnetic 
sites, an intermediate state with a total spin larger than 
that of the ground state has to be generated at one of the 
magnetic ions. Thus, there has to exist at least one spin- 
free orbital (empty or doubly occupied) at this ion, which 
is the case for most 3d ions except those in high spin d 5 



this expression can be written as 



T FE,1 _ jFE,X 
J ClC2 ~ <7 ClC2,r, 2 



J, 



FE,1 



ClC2,fl 



= FE 1 [(t fl , 2 ) 2 + (t C2 



(13) 



J FE contains two different hopping contributions, each of 
which is independent of the orbital £ of that site where 
the Hund interaction takes place. J FE is non-positive 
and leads to a parallel orientation of the spins Si and S2. 
This mechanism has been proposed by GKA to explain 
the occurrence of ferromagnetic coupling in non-metallic 
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compounds.^ Yet, since \FEi\ it is necessary that 

ti^ of Eq. H disappears or is very small at least. 

For case FD we consider a system of two equivalent 
magnetic ions Mi and M% with one singly occupied or- 
bital £ and one doubly occupied orbital /i each, see Fig. [j] 
c). As above, the ligand is empty. 

In contrast to case FE there are now 22 (= 2 • 11) 
channels, instead of two, all leading to a ferromagnetic 
coupling. The main difference between the channels is 
whether the spin of £ or one of the spins of fx first hops 



to the ligand. Putting 

J^' 1 = Jh 
C1C2 " 



A^ and A p 



1 



A c A M (A c + A„) 

I (± 

(A^(A C + A M ) VA M 

1 

(A M )2(C/ + A C -A M 



1 



V A c + A^ 
1 



A c + A M 
2 1 



(14) 



(A C )2(A C + A^ 



„ U + A c - A„ 

M2J 2 + (*C2 *Ml) 2 



A further simplification is obtained when we put A 



A L . Then 



C 



C1C2 



1 



3 2 

a7 + u 

[(*& + (*Ca W 

tFD,1 , tFD.1 

J ClC2^2 + J ClC2,Pl 



(15) 



= FD l [(t Clt 
Here, we have put 



u \u 



(k, 



(*C2 



2 

f Al- 

At 2 



3U 



(16a) 



(16b) 



jFD,i j g non _p OS itive and therefore leads to a parallel ori- 
entation of the spins. The expression is quadratic in the 
effective hopping-integrals t and the factor FD% is pos- 
itive (U,Al > 0). The coupling constants J^^ 1 ^ are 
independent of the orbital £ of the Hund interaction site. 
We note that FDi > FE%. As a consequence, the ferro- 
magnetic interaction should be more favorable for cases 
of a doubly occupied d orbital than for those of empty 
ones. This has already been pointed out by Andersonfl 



B. super-exchange via TV ligands 

In the following we consider the interaction via TV un- 
occupied ligands with one orbital each. We will assume 
hopping only between nearest neighbors along the chain 



of TV + 2 atoms and we investigate again the cases A, FE 
and FD. 

Case A requires a perturbation expansion up to order 

2(TV+ 1). There are f 2 ^^ 1 ^ possibilities (channels) to 

exchange the spins. To take account of all channels we 
use a simple method to generate them. For this purpose 
each channel will be written as a sequence of numbers 



(pip 2 • .-Pn+i\p 



N+2 ■ 



■P2N+1 P2N+2) 



(17) 



which represents the order of the hopping. Here p\ stands 
for the hopping from M\ to the first ligand, P2 for the 
hopping from the first to the second ligand, pn+i for the 
hopping from the last ligand to M2, pn+2 for the hop- 
ping back from Mi to the TVth ligand, see Fig. ||. The 
arrangement of the numbers in this sequence has to fulfill 
the restriction that the left as well as the right half must 
contain numbers of increasing magnitude. In the case of 
an expression of sixth order the numbers 1 to 6 have to 
be distributed. Therefore a possible sequence is 



(124|356). 



(18) 



First a spin from M\ hops to the first ligand (1) and 
from there to li (2). Then the spin from Mi hops to li 
(3) such that the two spins meet at li. Next M 2 becomes 
singly occupied again (4) . The last two processes are the 
hopping from li to li (5) and then back to Mi (6). Not 
possible is the sequence 



(3 5 6|2 1 4) , 



(19) 



since the ligands are empty and thus the hopping from 
li to l\ cannot occur first. 

With the simplification A^ = Aq 2 we find 



Here 



T A,N _ a (TN \2 



(W - m 



(N-l), n 2 
Hklk+x) 



N 

fc=l 



(20) 



(21a) 



(21b) 



The structure of Eq. ^T] is the same as that of Eq. [)] for 
ligand. Again, J A ' N is the product of a non-negative 
hopping-term and a positive energy factor An. Thus we 
have anti-ferromagnetic coupling. Comparing An with 
A\ indicates that J A,N is reduced as compared to J AA , 
because t N << t 1 . For large TV this is partly compen- 
sated by An > A\. 

Case FE involving an empty M orbital, can be han- 
dled in the same way. Nevertheless the structure of the 
sequences is somewhat different since the Hund interac- 
tion has to be taken into account. The sequences contain 
2 TV + 3 numbers which are ordered as follows 
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{pi Pi--- Pn+i \pn+2 \Pn+s ■ ■ ■ P2N+3) ■ (22) 

The central number indicates when the Hund interac- 
tion at a specific magnetic ion (either Mi or Ah), takes 
place, see Fig. |^. There is an additional restriction for 
the ordering of the p\. Both the series p\ . . -Pn+i and 
Pn+3 ■ ■ ■ P2N+3 as well as the triple pn+i \pn+2 \pn+3 have 
to contain numbers of increasing size. For this reason 
there is only one channel left. When two ligands are 
present this is 

(123|4|567). (23) 
Then, the interaction constant is given by 

T FE,N _ T FE,N , jFE.N _ pp r/rw n2 , ,7N \21 

(24) 

With A^ k A n we have 

FE N — FEi — (25a) 



The structure of FDn can be made plausible by com- 
parison with An. The perturbation series that was 
needed to calculate J A has the order 2N + 2, when TV 
ligands are present, i.e. it contains products of 2N + 2 
factors. The corresponding ferromagnetic interaction ad- 
ditionally includes the Hund interaction in first order. 
This factor can be introduced at 2N + 1 positions into 
the product of (2iV + 2) factors. In case of a {/-channel in 
An, i.e. when the spins meet at a magnetic ion, there is 
only one term of order U~ 2 , i.e. for pn+2 — 2N+1-. The 
other 2N cases lead to terms of order (U A) -1 . Then 
there are N ligand channels, here each of the 2N + 1 
terms give the same contribution of order A~ 2 . 

We have restricted the discussion of the coupling via 
more than one ligand to the case of ligands of the same 
type, i.e. they are either all anions or all cations. The 
generalization to the case of mixed-type paths remains 
work to be done. 



C. super-exchange via various paths 



IN 



tl k l k + 1 tl N Vj 



(25b) 



fc=i 



"^CiC2 m 1S non "P os itive and, because of FEn = FE\, it is 



much smaller than jf^' 1 ^ for large N. since t N « t 1 . 
There is no compensation from an increasing number of 
channels. 

Case FD: Here we find a rapidly growing number of 
channels. In sixth order these are 142. For this rea- 
son we restrict our investigation to the case of equivalent 
(Ai k — A Vfe) ligands in the chain. 

The sequences introduced in Eq. |2^ can also be used 
here. Yet the restrictions are modified. The series 
Pi . . - Pn+i and pn+3 ■ ■ -P2N+3 have to contain num- 
bers of increasing magnitude, while the center triple 
Pn+i\Pn+2\pn+3 now contains numbers of descending 
size. 

In this situation the coupling constant is (using Aq ss 



A) 



T FD,N _ jFD.N 

J ClC2 ~ J ClC2,M2 



FD.N 



FD N [(t£ 



ClC2,Ml 

TN V2 , rlN 

H2> + t- r C2M 



(26) 



J 2 



Here 



FD 



N 



J_H_ 

U 



}. + 2N± + N(2N + l)^ 



N 



(N-l) 



N P-j 



A N 



(27a) 



(27b) 



jFD,n j g n on-positive as is J™' 1 . As it contains terms 
of order N 2 it will compensate the decreasing magnitude 
of the hopping terms t N , even more than for the case A. 



So far we have discussed the super-exchange via one 
path. But clearly several paths of the same order can ex- 
ist along which the spins can be exchanged. A schematic 
view of this situation is shown in Fig. [3| 

We distinguish between equivalent and non-equivalent 
paths. Paths p and p' are called equivalent when Q p = 
Q p ' holds, with Q p — A P N , FE P N or FD P N . The additional 
paths can arise by introducing further ligand orbitals or 
by adding further ligand atoms. 

The interaction constant for a multi-path system is cal- 
culated by combining the contributions of the individual 
paths. The exchange takes place in two different ways. 
The first one is the use of only one path, forth and back, 
as before. The second one is a. ring exchange using two 
different paths forth and backB 

Again we consider the direct coupling of two half filled 
orbitals for the cases A, FE and FD. Altogether we find 
the following expressions 



p 

E 

p,p'=i 



(^ 



N 



Jy ' Ci C2 Ci C2 



(28a) 



(28b) 



p,p'=i 

7N,p 7N,p' , 7N,p fN,p' 

kim L CiV2 ^ ''Cam Cam 



■£SV) = 5 E ( FD N + FD p N ) 



(28c) 



p,p'=i 

7N,p 7N,p' 7N,p 7N,p' 
r ClMa C ClM2 ~T C C2Mi I C2Mi 
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P denotes the number of paths. 

If Q p = Q Vp the corresponding expressions for Jq are 
quadratic in the effective hopping- integrals. 




FE 



N 



E 

\p=i 



7N,p 
l CiV2 



\p=l / 



(29a) 



(29b) 



FE 



C1C2 



V1C2 



C1V2 



^E^SU+^E^ 



FD 



M1C2 



C1M2 



The factors , denoting the number of singly occupied 
orbitals at M t , result from the assumption of an orbital- 
independent on-site exchange. 

We will use equation ( |32j ) to calculate the coupling 
constants of two magnetic ions with arbitrary (integral) 
occupation which are coupled via unoccupied ligands in 
the frame of a Hubbard-model. We note that we are re- 
stricted to the case that all contributing paths are of the 
same order. 



jgf(P) = FD 




N 



(29c) 

In all above cases the existence of several equivalent paths 
only leads to a new effective hopping T^ N = J2p=i ^af- 
In particular the sign of the coupling constant remains 
unchanged, though the effective hopping T and thus the 
coupling may vanish when an appropriate set of param- 
eters is chosen. 

In case where the different paths are inequivalent this 
can lead to a ferromagnetic coupling resulting from di- 
rect supjer-exchange as has been discussed in a previous 
paper 113 



D. total interaction 

The total interaction J of two magnetic ions with total 
spins Si and S2 is obtained from Eq. (||) as the sum of 
the interaction of the half-filled orbitals 



J - E J d?2' 
C1C2 



where 



J CiC2 - j ClC 2 + E J C^2,vi + E J Ci 

+ E J o<2^i + E J <SL 



FE 



(30) 



(31) 



,M2 ' 



Here runs over all relevant unoccupied orbitals of 
Mi, J2u t over t ne doubly occupied respectively. Then it 
follows 



C1C2 



T A + V J FE 4- V T R 



FE 

(2,m 



(32) 



Hi 



E J £?2^i + E J a 



FD 

C2,M2 



III. FERROMAGNETIC COMPOUNDS 

In this section we apply our theory to selected insulat- 
ing ferromagnets, i.e. to FU^CrCU and the Cr tri-halides 
CrCl3, CrBr 3 and Crl 3 . In none of these cases an obvi- 
ous argument along the GKA rules can be found for the 
occurrence of ferromagnetism. 



A. Rb 2 CrCl 4 

This compound is a quasi- two-dimensional, ferromag- 
net with a Curie temperature Tc = 52.4 K.liil Values pf, 

J11 = -1.04 10" 2 eV and 4^ = 3 10~ 4 are reported.B 

11 jy 

The moments lie in the CrCl2-planes and are almost 
paralleLBO 

The crystal structure is analogous to that of K2CUF4. 
The parent structure is of K2NiF4-type, a body centered 
tctcagonal lattice with a = b = 5.086 A and c = 15.72 
A£3 The Cl _ ligands form an almost perfect octahedron 
around the Cr 2+ ions; yet, because of the d A configura- 
tion, a strong Jahn- Teller distortion takes place in order 



to remove the degeneracy of the e g -orbitals, singly oc- 
cupied in Cr 2+ . The distortion leads to elongation of 
the ligand cage along the y and x axes, respectively, 
and correspondinghz_.to contractions in the zx and zy 
planes, respectively.E3 The same distortions are found in 
K2CUF4.EJ Yet there is one important difference. While 
for K2CUF4 (dp configuration of Cu 2+ ) with its occu- 
pied d holes of alternatingly z 2 — x 2 and z 2 — y 2 type, 
the GKA rules can be applied in an obvious manner (the 
hopping between those orbitals is identically zero by sym- 
metry, and ferromagnetism results), no such rules apply 
for the occupied e 9 -orbitals of the Cr 2+ ions which are 
of 3y 2 — r 2 and 3a; 2 — r 2 type, respectively. This orbital 
order has been established by neutron measurements. lifl 
Consequently, er-type hopping is possible between neigh- 
boring orbitals. In addition, there exists 7r-type hop- 
ping between the singly occupied i2 S -orbitals. All these 
hopping processes favor an anti-ferromagnetic coupling. 
Nevertheless, the planar coupling is ferromagnetic. We 
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note that K2CUF4 and Rb2CrCU also show ferromagnetic 
coupling along the c-direction. 

In the hole picture for Cr 2+ , one of the e 9 -orbitals is 
doubly occupied. The other one is singly occupied (as 
are the t2g-orbitals) and alternates due to the orbital or- 
dering. When we define the hole state by 



g{8) = cos9\x 2 - y 2 > +sin6»|3z 2 



r 2 >, 



(33) 



we have values of 9i >2 = ±f ■ Then we find for the anti- 
ferromagnetic and the ferromagnetic contributions to the 
intra-planar coupling 



C1C2 

- 1 -A x (Cl) (2 {pdn f + \ (pda) 4 



A 2 
^ci 
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Here, Q denotes the singly occupied d-orbitals xy, yz, zx 
and .9(6*1.2). For simplicity we have put = A, i.e. we 
have ignored the level splitting of the t 2g and e g orbitals. 
The total coupling is given by 



J11 = 



1 

A 2 " 

^ci 

-(- — 

8\U + A C i 

3 J H ( 1 '2 
2~U 



(35) 
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The first two terms represent J A , where the first one, 
containing the factor (pdir) 4 , represents the coupling due 
to the t 2g orbitals. When we assume the 'canonical' ratio 
pdn/pda » —0.4, the t 2g orbitals contribute less than 1/3 
to J A . Note that any mixed terms (pda) 2 (pdir) 2 etc. do 
not exist, since the e g and t 2g orbitals are decoupled by 
symmetry. The third term represents J F , which contains 
a factor Jh/U « —0.1. Yet, on the other hand, the ge- 
ometry of the a- hopping provides a factor 3/2 (compared 
to 1 /8 in term 2) . In addition, the large number of chan- 
nels in FD 2 leads to a reduction of the effective energy 
denominator by approximately a factor of 2. 

As a consequence, a net ferromagnetic coupling is ob- 
tained. This result is rather stable, when energy and 
hopping parameters are varied within reasonable limits. 
We note that this situation may change, when further 
octahedral distortions occur which mix the t 2g and e g 
orbitals. Probably, there exists an easy axis single ion 
anisotropy, favoring spin alignments either along the y 
or the x axis. The ferromagnetic nearest neighbor cou- 
pling will then produce the almost complete alignment of 
spins. 



For the related case of KCrF3, with the same struc- 
ture of Cr planes as in PJ^CrCU suggesting the same 
orbital structure of the Cr ions, the puzzle of ferromag- 
netic intra-planar coupling has been studied previously 
by other researchers. t3't3 They have assumed direct hop- 
ping between the magnetic sites and have suggested that 
values of 9 m j should occur, in contrast to the findings 
in Rb 2 CrCl 4 , which yield 6 = f .0 Note that the neglect 
of ligand orbitals strongly reduces the number of ferro- 
magnetic channels - there is only one left as FDq = FEq 
holds. The coupling J A , which vanishes only for 6* = -|, 
is relatively large for 9 — -| and wins over FD . A fur- 
ther, even more questionable consequence is the fact that, 
neglecting the ligands, the energy change due to the oc- 
tahedral distortion is of the same order as the magnetic 
interaction. When the ligands are included, the energy 
gain due to the Jahn- Teller distortion dominates over any 
magnetic energies. An analogous d 4 system with planar 
ferromagnetic coupling is LaMn03 . Again we think that 
simplified coupling estimates neglecting the ligands are 
not appropriate. E3 

An analysis of the inter-planar coupling has not been 
performed. On the one hand our formalism is not ap- 
plicable to the situation with occupied and unoccupied 
ligands. On the other hand the number of spins is large 
(16 in the hole picture) such that a numerical calcula- 
tion is very extensive. But we expect a similar behavior 
as in K2CUF4, where a reduction of J A due to destruc- 
tive interference with a possible sigii|-ehange and a more 
or less unmodified J F was obtained.113 Furthermore, for 
Rb 2 CrCi4 the spin- free orbital is doubly occupied and 
therefore the case FD holds, which always leads to larger 
ferromagnetic coupling than the case FE. 



B. CrHa 3 (Ha = Cl, Br, I) 



The Cr tri-halides with rhombohedral symmetry (R3) 
contain two formula units per trigonal unit cell. The 
two Cr atoms occupy equivalent positions, so do all six 
halide atoms. They form triple layers Ha-Cr-Ha which 
are stacked along the c-axis. In the central Cr layer, the 
Cr atoms form a honeycomb lattice like a single layer 
of graphite. It is helpful to view this plane as a closed 
packed triangular arrangement of two Cr atoms (I and 
II) and of a vacancy V in the center of the honeycomb. 
Both the Cr atoms and the vacancy are surrounded by 
edge sharing octahedra of halide atoms. As a result, the 
two halide planes above and below are closed packed tri- 
angular layers with a stacking sequence A-B (see Fig. ||) . 
The triple layers are now repeated along c in such a way 
that the vacancies V successively shift by |(a — b + c) , 
where a and b are the basal plane lattice vectors (a = b, 
a-b = — i) and c is the c axis vector of the corresponding 
hexagonal unit cell. The motion of V from positions A 
to B to C to A (see Figs. |] and ||) leads to a packing of 
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the halide planes along c of the form A-(Cr)-B-(empty)- 
A-(Cr)-B-(empty)-A-(Cr)-B etc. As a consequence, the 
interlayer space B- (empty)- A also consists of halide octa- 
hedra. One can define a perfect octahedral closed packed 
halide lattice by requiring that all octahedra are regular 
(rhombohcdral angle a — 60°). This would lead to the 
same separation of A and B planes whether or not a Cr 
layer is intercalated. It is clear however that the inter- 
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planar distances are quite different, see Tab. III. The lat- 
tice parameters of the Cr tri-halides are given in Tab, p, 
some relevant atomic separations are listed in Tab. Ell. 



Note that for our calculations of magnetic coupling con- 
stants it was advantageous to start from an 'ideal' rhom- 
bohcdral lattice by assuming regular octahedra around 
the Cr atoms. In the real lattice, there exists a compres- 
sion of the octahedra enclosing the Cr plane. As a result, 
the octahedra surrounding the empty layer are consid- 
erably elongated along the c direction. There are two 
further types of distortion allowed by the R3 symmetry, 
which, however, appear to be of minor effect. 

1) The Cr planes buckle, i.e. , the Cr(I) and Cr(II) atoms 
lie respectively above and below the plane, as defined by 
the vacancy positions (deviation from 1/3 of the lattice 
parameter u, see Tab. |ll]). 

2) There are some minor shifts in the halide positions 
(lattice parameters x deviating from 1/3 or 2/3 values 
and y deviating from zero). These shifts result partly 
from a relaxation of the octahedra around the vacancy 
V and partly from the Cr buckling. Note that a slight 
rotation of the V octahedra follows, which should lead to 
optical activity along the c axis. 

Of the three compounds, the lattice structure of CrBr3 
has been investigated most extensivelyJlalij Here, the in- 
ternal lattice parameters x,y,z,u are known, while for &I3 
these quantities have not been determined. For CrCl3, 
Morosin and Narath report a low temperature transition 
from an orthorhombic to the R3 near T = 225 K and give 
values for x,y,z,u.E9 The difference of the two structures 
are different stacking arrangements of the Cr tri-halide 
planes. To our knowledge, the CrCl3 structural phase 
transition has not been confirmed by other groups. There 
is evidence that CrCl3 remains orthorhombic_to low tem- 
perature into the anti-ferromagnetic phase.EU 

All three compounds exhibit strong intra-planar ferro- 
magnetic coupling. This feature is to be expected from 
the GKA rules, as the Cr-Ha-Cr angle for the nearest 
neighbor (fourth order) super-exchange paths (leading to 
the coupling constant J\ ) is rather close to 90° - it is 94° 
for CrCl 3 and 93° for CrBr 3 . 

The inter-planar coupling is anti-ferromagnetic for 
CrCl3, while it is ferromagnetic for the other two ma- 

16.8 K, 



terials. 
1 c — 



-<Cl 
■-N 



The critical temperaturns_£ 
32.5 K and T£ = 68 K.E3 Note that an easy 
plane anisotropy is reported for CrCU. while easy axis 
anisotropies occur in CrBr3 and Crhj.Ej 

From analyses of magnetic susceptibility data intra- 
planar coupling constants Jj| have been obtained. The 



-4.1 meV, Jf r 



values reported in Ref. [| are J,j 

—6.4 meV and J,| = —10.5 meV. On the other hand, 
Samuelsen et al. have carried out neutron measurements 
of spin wave dispersion curves for CrBr3.Ha They have 
been able to extract several near neighbor coupling con- 
stants, in particular for the inter-planar coupling. For 
the intra-planar coupling, they report J\ = —7.6 meV 
and J3 = 0.13 meV (next nearest neighbor, sixth order 
super-exchange). Note that \J^ r \ < \J\+2J- i \. Since 



the J 1 value of Ref. 18 probably is the most reliable ex- 



perimental quantity available, we have used the relation 
J\ = 1.2JP r to scale Jp l and J| accordingly and have 

assumed Jj 7 ' = -4.8 meV and J{ = -12.4 meV. Note 
that we use the definition of J of Eq. |5|, leading to a multi- 
plication by a factor of nine of the values originally given 
by Refs. | and 

There exist three different inter-planar coupling con- 
stants of sixth order super-exchange, labeled J 2 , J4 and 
J5 by Samuelsen et al. J2 is the vertical Cr(I)-Cr(II) cou- 
pling constant, J4 represents another I-II coupling, while 
J5 stands for I-I or II-II couplings (see Fig. |). The val- 
ues given by Ref. [l^ are J2 = 0.22 meV, J4 = —0.06 meV 
and J 5 = -0.12 meV. 

In the following we present results of calculations for 
the coupling constants J\, J2, J4 and J5 using the the- 
ory of Section II. We have not included J3, as we mainly 
focussed on the question of inter-planar coupling. 

The electronic structure of the Cr ions depends on the 
ligand position. Because of almost perfect octahedra the 
three singly occupied states are almost t2 g -like and the 
two doubly occupied ones have mostly e g character, as 
can be found from diagonalizing the d-states in the ligand 
field. Note that Cr 3+ is not a Jahn- Teller ion. The three 
empty p-states of the halides are treated as degenerate. 

The intra-planar coupling of two nearest Cr ions takes 
place via two ligand ions, see Fig. ||. From Eq. ( p2| ) the 
exchange interaction is given by 



■h = Ai(Ha) ^(T C \ C2 ) 2 + 6FDi(Ha) 

C1C2 



Cm 



(36) 



and and the energy 



VI. The effective hop- 



with the effective hoppings 
factors Ai and FDi, Tabs 
pings are calculated corresponding to their definition in 
Eq. ( p9| ) . For the energy denominators, we have neglected 
the energy difference between e g and igo like states. 



The hopping parameters from Tab. IV are found by a 
fitting process, where the charge-transfer energy A was 
determined by adjusting J^ a to the experimental values 
as given above, while U and Jh were kept fixed. We note 
that the order of magnitude and tendency of the A-p a and 
the value for U agree with those given by Zaanen.El These 
values are also consistent with those found for Rb2CrCl4. 

The calculation of the inter-planar coupling is very 
complicated. On the one hand up to six ligands with 
three orbitals each are involved in lowest order exchange. 
On the other hand there are three occupied ground state 
orbitals. 
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For the calculation of J 2 , J a and J5 we need the di- 
rect (here anti-ferromagnetic) and the indirect (ferromag- 
netic) coupling constant in sixth order. Each of these 
couplings is mediated via several paths. Nevertheless, be- 
cause we have assumed that all ligand orbitals are equiv- 
alent, this only leads to the introduction of an effective 
coupling as was discussed above. 

The vertical inter-planar coupling J 2 is mediated via 
six halide ions with three orbitals each. The halide ions 
form an (empty) octahedron elongated along the c axes, 
see Fig. ||. Three halide ions are ligands of the lower 
Cr ion and the other three those of the upper one. The 
couplings J4(I-II) and J5 (1-1,11-11) are mediated via four 
halide ions where to each Cr belong two ligands, see 
Fig. § 

For the coupling constants we get from equations (|2E 
and @ 



J = A 2 (Ha) J2KcJ 2 + 6FD 2 (Ha)J2(T?») 2 - 



(37) 



C1C2 



Thus we have to determine the effective hopping integrals 

T (i(2 ' T (i»2 and T ^i<2 and tne ener gy factors A 2 (Ha) and 
FD2{Ha) (independent of the path). 

The terms T^p depend on the individual Cr-halide and 
halide-halide hopping integrals. In particular, the ratio 
(ppa) / (ppir) of the halide-halide hopping integrals may 
modify the relative magnitudes of the J 2 , J a and J5 cou- 
plings. The values for Tj^ for the inter-planar couplings 
are shown in Tab. [v| and the resulting energy factors as 
well as the values of the J 2 , J a, J5 coupling constants 
are given in Tab. [VT| . For CrBr3 these numbers compare 
well with the neutron data of Ref. |l8], though the ratio 
J4/J5 is not very good. Note that the 'vertical' J 2 always 
is relatively large and anti-ferromagnetic. In the case of 
CrBi - 3 andn Crl3, this coupling is over-compensated by 
the ferromagnetic J a and J5 couplings. This is also found 
for CrCbj unless the ratio \ (ppa) / (ppn)\ was shifted from 
the usual range of 3 - 5 to an extremely large value of 



36, such that J c = J2i= 



2 4 5 z i Ji > 0. If so, the inter- 
planar coupling is probably anti-ferromagnetic for CrCl3 
with this set of parameters. (But note, that J c can not be 
identified with the inter-planat-icoupling in a model 
with two couplings J11 and Jl£3) Nevertheless, we think 
that \{ppa) / {ppir)\ = 36 is outside the range of physically 
meaningful parameters. Thus we suspect that there may 
indeed exist a difference in the low temperature struc- 
tures of CrCl3 and the other two materials. We finally 
note that the dominant magnetic coupling constant of 
these materials is the intra-planar J\ (4th order super- 
exchange). As we pointed out in Sec. II A, the -FlD-type 
coupling should cause a minority spin polarization on 
the halide ligands, as observed in the neutron study of 
spin density distributions in CrBr3 by Radhakrishna and 
Brown J13 



IV. CONCLUSION 

Nonmetallic transition-metal oxides and halides ap- 
pear to be the most prominent members of the family 
of Mott-Hubbard insulators. In this paper we have gen- 
eralized the expansion of extended Hubbard models by a 
Rayleigh-Schrodinger perturbation series into Heisenberg 
models, with the goal to determine general expressions 
for the super-exchange coupling constants. Each part of 
this series that corresponds to an exchange of the spins in 
the half-filled orbitals (1 and £2 of two magnetic ions Mi 
and M 2 gives a contribution to the (isotropic) coupling- 
constant J = X<iC2 °^CiC2 °f the Heisenberg Hamiltonian 



H = J 



Si ■ S2 

4(Si)(S 2 )- 



(5) 



Here Si denotes the total spin at the magnetic site Mi. 
This form of the effective spin Hamiltonian follows from 
the assumption of a relatively strong Hund interaction at 
the magnetic sites such that Si = S^ . 

The couplings Jc,iQ 2 are built from two different parts. 
Either the exchange takes place involving only the d or- 
bitals Ci and £2- This 'direct super-exchange' interaction 
is anti-ferromagnetic (A) in general. Or 'indirect super- 
exchange' occurs, which involves an additional, empty 
(E) or doubly occupied (D), d orbital. In this case, 
because of the Hund interaction, the coupling is nega- 
tive in general. These general results correspond to the 
Goodenough-Kanamori- Anderson (GKA) rules as stated 
in the introduction, yet quantitative results - and this 
may include the sign of the coupling constant - depend 
strongly on details of the involved exchange process. The 
simplest cases are those when only one specific path 
via ligands occurs. Then, each hopping integral enters 
quadratically the expressions for the coupling constants 
J. Therefore, the sign of J is given by kind of super- 
exchange, direct (+) or indirect (— ). The magnitude 
of J depends on i) an effective energy denominator Eq 
and ii) on an effective term t, such that J® = P/Eq. 
Here Q — A, FE, FD represents the different possibili- 
ties of super-exchange and i — (Yl^^i tk)/A N , where N 
denotes the number of ligands involved in the exchange 
and the tk the corresponding hopping integrals. A is a 
characteristic energy of an excited state with tk/ A <C 1 . 
Thus t decreases rapidly with increasing path length N . 
1/Eq, which we also call Qn, is strongly influenced by 
the number of exchange channels. This number depends 
on the specific super- exchange A, FE, or FD, each of 
which exhibits a different dependence on TV. FEn is in- 
dependent of N, there are always only two channels. An 
is of the form A N = 2 (l/U + N/A) and FD N includes 
in addition a term quadratic in N. For N = 1 one finds 
FEi < FD\ < Ai, mainly because FEn and FDn con- 
tain an additional factor Jh/U. With growing number 
of ligands we may find cases An < FDn- As a conse- 
quence, when N > 1, the interaction of ions with more 
than half-filled d-shells should exhibit a bigger tendency 







towards ferromagnetic coupling than those with less than 
half- filled d-shells (hole picture, empty ligands). 

The exchange can also take place via several (P) in- 
terfering paths. The structure of a coupling constant 
resulting from P paths with N ligands each is J®(P) = 

h J2pp'=i (Qn+Qn)^ t p ■ Here p and p' indicate the var- 
ious paths. The energy denominators are the same 
as in the case of P = 1, yet the hopping factors no longer 
enter quadratically into the expressions for the effective 
hopping integrals. In the presence of inequivalent paths 

(Q P N ^ Q P N ) this leads to the possibility of a ferromag- 
netic coupling of two spins by direct super-exchange as 
was discussed in Ref. [lC]. When all P paths are equiva- 
lent, the J Q are again quadratic in the effective hopping, 
i.e. the sign is fixed, but the magnitude may be strongly 
influenced by interference. This behavior always occurs 
for case FE, since FE^ is independent of N. 

From the above considerations we see that even the 
sign of the coupling cannot be guessed except in the sim- 
plest cases, even if the orbital order at the magnetic sites 
is known. For better quantitative estimates, a reasonable 
knowledge of hopping integrals, of charge transfer gaps 
and of U parameters is required. These numbers are ac- 
cessible, in principle, via total energy calculations using 
constrained density functional theory, see e.g. Ref. |2J. 
Our present studies are based on such evaluations for 
La 2 Cu0 4 . 

In order to test the quality of our treatment, we have 
investigated some ferromagnetic Cr compound insula- 
tors, where the GKA rules cannot be applied in any 
simple fashion. In RbaCrCLi, with the same structure 
as K 2 CuF4, orbital order of the Cr 2+ ions occurs be- 
cause of the Jahn- Teller effect. However, for the first 
nearest neighbor (INN) intra-planar coupling, the direct 
super-exchange J a does not vanish, in contrast to the 
case of K2CUF4. Nevertheless, we find \Jp D \ > J\, for 
various reasons as discussed above and in detail in chap- 
ter 3. Also the quantitative result agrees well with ex- 
periment. Concerning the inter-planar coupling, we be- 
lieve that Rb2CrCU is rather, similar to K2CUF4 which 
has been studied previously.E9 We note that some of the 
N = 2 super-exchange paths involve both empty (CI - ) 
and doubly occupied (Rb + ) ligands. Such cases, which 
may occur in ternary compounds, are not covered by our 
present theoretical treatment. 

The Cr 3+ ions of the Cr tri-halides CrCl 3 , CrBr 3 and 
Crl 3 are not Jahn- Teller ions and, thus, do not exhibit 
orbital ordering. The INN intra-planar coupling involves 
two equivalent N — 1 paths with Cr-Ha-Cr angles close 
to 90°. The FD coupling dominates the A coupling, both 
due to the 90° angles and, also because of the very many 
FD channels. The inter-planar coupling takes place via 
three different N = 2 contributions, J2, J a and J5. Each 
of the couplings is mediated via a variety of paths. Using 
reasonable hopping parameters we get good agreement 
with the values J2, J4, J5 as measured by neutron scat- 
tering in the case of CrBr 3 . Further we obtain ferromag- 



netic coupling for both CrBr 3 and Crl 3 . In the case of 
CrCl 3 we suspect that the low temperature structure is 
not R3. 

We come to the conclusion that the perturbational cal- 
culation of the coupling constants in the non-metallic 
3d-transition-metal compounds is an appropriate instru- 
ment, qualitatively as well as quantitatively. The quan- 
titative agreement may appear somewhat surprising. It 
has been pointed out earlier that the higher order con- 
tributions of the perturbation series can have a similar 
magnitude as the lowest one.E3 Under such circumstances 
we would expect that the coupling between next nearest 
neighbors becomes important. Ex-perimentally, however, 
ratios J2/J1 ~ 1CU 2 are reported£l An additional conse- 
quence should be that terms of the form (Sj • Sj)(Sfc • S;) 
can not be neglected. We are not aware of any experi- 
mental evidence for the importance of such 4-spin terms. 
The question remains how the good agreement between 
experimental value and the lowest order results can be 
understood. 



FIG. 1. Spin interaction by super-exchange via one empty 
ligand with one orbital, a) anti-ferromagnetic case, b) and 
c) the ferromagnetic cases. The double arrow represents the 
direct exchange Jh that leads to the second Hund's rule, b) 
and c) are differ by the occupation of the spin-free orbital at 
the magnetic ion. It is not possible to transform one into the 
other by a particle-hole transformation. 
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FIG. 2. a) Direct interaction on a linear chain between two 
half-filled orbitals £1 and (2 at the sites Mi and M2 via N un- 
occupied ligands with one orbital each. Only nearest neighbor 
hopping is considered, b) Indirect interaction between two 
half-filled orbitals £1 and £2 at the sites Mi and M2 involving 
one unoccupied orbital i) via N unoccupied ligands with one 
orbital each which build a linear chain. Only the hopping 
between nearest neighbors is allowed. 
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FIG. 3. Super-exchange via three different paths with three 
ligands each. 



FIG. 5. Crystal structure of the tri-halides CrCLj, CrBr3 
and CM3. The hexagonal unit cell contains four Cr-planes. 
Large circles denote the halide ions the smaller circles the Cr 
atoms. Clusters corresponding to the couplings Ji, J2 and 
J4, see Fig. [| are indicated. Also shown are the position of 
a ligand octahedron and the corresponding local coordinate 
system of a Cr ion. The distances d, 1 and s are given in 
Tab. 



Ill 
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FIG. 4. Projection of a Cr-plane with the corresponding 
ligands on the a6-plane. The small circles denote the Cr-ions, 
the large ones the halides. The darker ligand ions lie below the 
Cr-plane the lighter ones above. The Cr planes are stacked 
such that the Cr hexagon center succesively lies at points A, 
B and C. 




FIG. 6. Near neighbor coupling constants Ji, J2, Ji and 
J5 in CrHa3 . There exist three couplings Ji , one J2 , three J4 
and six J5 per Cr. 



TABLE I. Electronic parameters for RL^CrCU (in eV). 



(pda) (pdn) U Jh Acz 
-1.0 0.4 4 -0.52 4.15 



TABLE II. Lattice parameters for CrCl 3 (Ref. |l|), for 
CrBr 3 (Ref. [jj) and for Crl 3 (Ref. §(]). For Crl 3 the x, y, 
z and u values are approximated by the 'ideal' values, 2/3, 0, 
Zi = ao/2\/3co , and 1/3, respectively. 



CrCla CrBr 3 Crl 3 

5.94A 6.27A 6.86A 

17.33A 18.21A 19.88A 

0.6507 0.6523 0.6666 

-0.0075 0.0012 0.0 

0.0757 0.0786 0.0813 

0.3323 0.3339 0.3333 



a 

Co 
x 

y 

z 
u 
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TABLE III. Relevant atomic separations for the CrHa3 
compounds (in A, for definition, see Fig. ^). Because of the 
distortions u 7^ 1/3, y 7^ etc. there exist two different 
distances d (Cr-Ha) and I (Ha-Ha). Also given are the sepa- 
rations of the Ha planes n (Cr enclosed) and r 2 (empty site 
enclosed) 





CrCl 3 


CrBr 3 


Crl 3 


s 


3.43 


3.62 


3.96 


<h 


2.37 


2.50 


2.80 


<h 


2.33 


2.50 


2.80 


h 


3.69 


3.80 


4.09 


h 


3.80 


3.87 


4.09 


ri 


2.62 


2.86 




'•2 


3.15 


3.21 





TABLE VI. CrHa3: Effective energy denominators A and 
FD (in eV _1 ) and coupling constants (in meV). 





CrCl 3 


CrBr 3 


Crl 3 




1.10 


1.18 


1.24 


u r U\ 


QzL 


111 


1 98 
— ± .^0 


A 2 


1.64 


1.77 


1.91 


6FD 2 


-2.28 


-2.76 


-3.36 


Ji 


-4.8 


-7.6 


-12.4 


■h 


0.059 


0.26 


0.31 


Ji 


-0.007 


-0.08 


-0.33 


J 5 


-0.006 


-0.09 


-0.33 



TABLE IV. Electronic parameters (in eV) for CrHa 3 . The 
hopping integrals are given for the 'ideal' structure. Dis- 
tance variations are included via exponentials e' ,3 " ro_r ^ with 
= 1.5 eV/A, (5l d = 2 eV/A, = 2 eV/A and = 2.5 
eV/A, where tq denotes the distance in the ideal structure. 





CrCl 3 


CrBr 3 


Crl 3 


A 


3.8 


3.3 


3.0 


U 


3.5 


3.5 


3.5 


Jh 

u 


-0.1 


-0.1 


-0.1 


(pda) 


-1.15 


-1.12 


-1.06 


(pdn) 


0.66 


0.69 


0.70 


(ppa) 


0.05 


0.10 


0.12 


(ppn) 


-0.0014 


-0.019 


-0.04 



TABLE V. CrHa3: Values of various effective hopping 
squares for the couplings J\, J 2 , J4 and J 5 (in 10~ 4 eV 2 ). 





CrCl 3 


CrBr 3 


Crl 3 






Ji 






9733 


7670 


4802 




12211 


8862 


5506 






J2 






84 


359 


444 




12 


98 


182 






Ji 




Scic 2 ( t2 ic 2 ) 2 


18 


79 


107 




18 


88 


140 






J5 






11 


69 


107 




13 


85 


140 
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